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Abstract Nuclear magnetic resonance (NMR) experi- 
ments are described for two granular-flow systems, the 
vibrofluidized bed and the gas-fluidized bed. Using pulsed 
field gradient, magnetic resonance imaging, and hyper- 
polarized gas NMR, detailed information is obtained for 
the density and motions of both grains and interstitial 
gas. For the vibrofluidized bed, the granular temper- 
ature profile is measured and compared with a first- 
principles formulation of granular hydrodynamics. For 
the gas-fluidized bed, dynamic correlations in the grain 
density are used to measure the bubble velocity and hy- 
perpolarized xenon gas NMR is used to measure the 
bubble-emulsion exchange rate. A goal of these measure- 
ments is to verify in earth gravity first-principles theories 
of granular flows, which then can be used to make con- 
crete predictions for granular flows in reduced gravity. 
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Fig. 1 Routes to reliable predictions of granular-media 
properties in reduced gravity. Discrete-element (DE) mod- 
eling (sometimes called molecular dynamics (MD) modeling) 
has the acceleration of gravity g as a direct input. Similarly, 
microscopic theories of granular media depend explicitly upon 
g. Both types of calculation should be tested by physical 
experiments in earth gravity, before more difficult reduced- 
gravity experiments are carried out. 
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1 Introduction 

Granular flow and fluidized-bed properties in reduced 
gravity and interstitial gas pressure will be important 
for processing, storage, and other applications during lu- 
nar and Martian missions. However, current engineering 
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practice for granular flows is based on correlations of 
Earth-gravity data and so cannot be extrapolated with 
any certainty to reduced gravity. An alternative path is 
to use experiments to validate granular flow equations 
derived from fundamental physical principles, in which 
the acceleration of gravity appears as an explicit input 
parameter. To the extent that this is possible, reasonable 
predictions for reduced gravity granular systems can be 
made, although it will remain necessary to test the pre- 
dictions in actual microgravity experiments (Fig. ^| . 

Nuclear magnetic resonance (NMR) enables nonin- 
vasive measurements to be made on dense, opaque sys- 
tems such as granular media. We are engaged in a pro- 
gram of experiments to build up detailed pictures of 
dense, three-dimensional granular systems in Earth grav- 
ity using NMR methods combining MRI, PFG-NMR, 
and hyperpolarized-gas NMR. Using these methods we 
are able to measure density profiles, and grain as well as 
gas motions over millisecond to second time scales. We 



Fig. 2 Typical NMR pulse sequence used for grain mea- 
surements in granular flows. The vertical lines indicate radio- 
frequency pulses, and the rectangles indicate gradient pulses 
of variable strength. Using this type of sequence, it is possi- 
ble to measure the joint distribution of grain motions during 
the interval At with pim resolution, and grain positions with 
sub-mm resolution. The entire sequence can be completed in 
as little as 2 ms, although it must be repeated many times to 
build up the distribution information. 

will describe our studies of two granular flow systems: the 
vibrofiuidizcd bed and the bubbling gas-fluidized bed. 



2 NMR techniques for granular flows 

Broadly, two requirements must be met to carry out 
NMR experiments on a flowing granular medium. First, 
a system must be set up in the desired granular flow 
state (e.g., bubbling gas-fluidized bed) which also gives 
a substantial NMR signal for the subsystem to be mea- 
sured (grains or gas). Second, suitable NMR protocols 
must be devised to extract the desired information such 
as density, velocity, diffusion, exchange between phases, 
etc. The ways in which these two requirements are met 
are quite different for grain and gas NMR, yet sometimes 
it is possible to carry out grain and gas measurements 
on nearly identical systems. 

2.1 Grain NMR techniques 

Due to the rapid molecular tumbling, fluids give NMR 
signals much more suitable for imaging-type studies than 
do solids . Botanical seeds, which typically contain liq- 
uid oil that gives a strong NMR signal, have been one 
popular type of sample for granular-media studies [SIT^|. 
To allow NMR studies over a wider range of grain sizes, 
small amounts of organic liquids can be adsorbed onto 
substrates such as porous catalyst particles [2l|. In the 
experiments described here, both types of grains were 
employed. 

For grain-phase studies, we typically use NMR pulse 
sequences that can encode both the motion of the grains 
and their positions (Fig. [2J . Using strong gradient pulses 
(up to 1000 G/cm), it is possible to resolve grain motion 
over times of one millisecond or less. In this type of study, 
the individual grains are not imaged. Rather, the joint 
probability distribution over all grains is measured as a 
function of grain position and grain displacement in a 




Fig. 3 Joint distribution p(Az, z) of vertical displacements 
Az in At = 1.38 ms and heights z, f or the grains in a vi- 
brofluidized bed of mustard seeds ll ill . Near the bottom of 
the bed (small z) the distribution is highly skewed due to 
collisions with the vibrating bottom wall, while at large z the 
distribution is nearly Guassian. 

specified time interval Q ■ Figure shows an example of 
this type of data. Note that for our studies of systems 
with continuous, rapid grain movement both the NMR 
techniques used and the type of data obtained are quite 
different from ealier studies of quasistatic flows Q. 



2.2 Gas NMR techniques 

Under equilibrium conditions, the NMR signal from gasses 
is very weak (due to their low density), often making 
gas-phase NMR studies difficult. By using hyperpolar- 
ized gas, which has been shown to be practical for 3 He 
and 129 Xe, the NMR signal can be increased by a factor 
of order 10 4 |2^. Figure 0] shows a simplified diagram 
of the apparatus we use to study the gas motion in a 
gas-fluidized granular bed. In this apparatus, xenon gas 
is polarized by spin exchange with rubidium vapor in a 
laser optical pumping cell. The hyperpolarized xenon is 
then used as the fluidizing gas for a small granular bed 
located within the NMR magnet. 



3 The vibrofluidized bed: Granular 
hydrodynamics 

The first granular system to which we have applied these 
techniques is the vibrofluidized bed. Mechanical vibra- 
tion is a fundamental and widely used technique for trans- 
forming granular media from solid-like to fluid-like states. 
In addition, vibration is often used as an aid for other flu- 
idization methods such as gas flow or gravity-driven flow. 
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Fig. 4 Simplified diagram of apparatus for using hyperpolarized xenon to study gas motion in a fiuidized granular bed. 



Questions one would like to be able to answer from ba- 
sic physical principles include: What amplitude and fre- 
quency of vibration are required for fluidization? What 
are average grain density and grain speeds as functions 
of height within the bed? 

One theoretical approach has been to treat the indi- 
vidual grains like molecules in a (pseudo) fluid, and to 
bring to bear methods of statistical physics such as ki- 
netic theory and hydrodynamics |j| El El UB U2] ■ This 
approach is most tractable (and more likely to be valid) 
in the "rapid granular flow" regime, in which the tran- 
sient two-body collisions between the grains dominate 
over longer- lived and many-particle interactions |3(. For 
a vibrofluidized bed, this regime occurs for strong, high 
frequency driving of relatively shallow beds . 

Several earlier studies had used video techniques to 
measure grain motion statistics in q uasi-two dimensional 
vibrofluized beds 0; El El El H3j- Some of these two- 
dimensional studies found strong deviations of the grain 
velocity distribution from the nearly-Boltzmann distri- 
bution that might be expected from the statistical mod- 
els. Using NMR, we have been able to noninvasivly probe 
the grain velocity distribution in a three-dimensional bed. 
For the conditions of our experiments we find that the 
grain velocity distribution is close to the Boltzmann dis- 
tribution (Fig. EJ, validating a statistical approach. Be- 
yond this, we are able to make a detailed joint fit of the 
measured density and granular temperature profiles, to a 
first-principles calculation Q of granular hydrodynamic 
coefficients (Fig. EJ). 

The hydrodynamic equations used to generate the 
theory curves in Fig. ^explicitly reference gravity (using 
a Bernoulli relation, see Ref. for details). After val- 
idating the hydrodynamic theory in Earth gravity by a 
series of fits like that shown in Fig.^J it is reasonable and 
easy to extend the theory to reduced gravity by recalcu- 
lating it with reduced g. The output of such a calculation 
could be a prediction for the density profile in a granu- 
lar bed vibrofluidized in lunar or Martian gravity. This 
is a (perhaps rather simplistic) example of the program 
outlined in Fig. ^ 
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Fig. 5 Distribution of horizontal grain displacements in a 
vibrofluidized bed under the same conditions as Fig. |3] The 
NMR observation time At was 1.38 ms, consi der ably shorter 
than the mean grain collision time « 6 ms |TJ|. Therefore, 
the measured displacement distribution should give a good 
approximation to the grain velocity distribution. The curves 
show fits to a Gaussian function with and without power-law 
stretching of the exponent. The distribution deviates only 
slightly from an unstretched Gaussian function, which would 
describe the distribution of velocities of molecules in a molec- 
ular fluid. 



4 The gas-fluidized bed: Grain and gas NMR 

Gas-fluidized granular beds are widely used in industry 
for catalytic cracking, efficient combustion, and many 
other large-scale processes El El It is reasonable to 
imagine that gas fluidization (or pneumatic transport of 
granular material, which is closely related) would play a 
prominent role in in-situ resource utilization on the moon 
or Mars, yet the behavior of a gas-fluidized bed in re- 
duced gravity is almost completely unknown at present. 

To promote mixing, gas-fluidized beds are typically 
used in the bubbling regime, in which large bubble-like 
voids devoid of grains rise through the bed at a velocity 
much larger than the mean gas velocity. From a theo- 
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Fig. 6 Density n* — nd 3 (top) and granular temperature 
T* = T/mgd (bottom) as functions of height z* = z/d, 
for a vibrofiuidized bed under the same conditions as Fig. 
(the starred variables have been nondimensionalized by the 
grain diameter and mass d, m and the acceleration of gravity 
g). The symbols show NMR measurements, while the curves 
show a fit to the granular- hydrodamics theory of Ref. 0- A 
notable feature (in both the experiment and the theory) is 
a "temperature inversion" — even though the energy flow is 
strictly upwards, the temperature falls with increasing height 
z* and then rises again at large heights. At heights above the 
temperature minimum, the flow of granular "heat" is from 
colder to hotter temperatures, an effect that is only possible 
in highly out-of-equilibrium systems. 



retical standpoint the bubbling gas-fluidized bed is ex- 
tremely complex, with significant motion over a wide 
range of length scales (from grain to bubble) and fluid 
dynamics with self-consistently moving boundary con- 
ditions. Microscopically- motivated theories of fluidized- 
bed behavior have been formulated (see, for example 
Ref. 0| ) yet to achieve closure these theories must make 
various unverified assumptions. For example, in a "two 
fluid" model it might be assumed that the grains and gas 
behave like two interpenetrating Newtonian fluids [12 ]. 
Yet, at present it is not established theoretically or ex- 
perimentally that a dense granular medium is ever well 
described as a Newtonian fluid. 

In common with other granular media, the gas-fiuidized 
bed is difficult to probe experimentally in a non-invasive 
manner. To better understand the dynamics of the bub- 
bles and their interaction with the dense ("emulsion") 
phase of the bed, many techniques have been tried in- 
cluding x-ray photography and tracer-gas injection 



Here we show preliminary data for two projects we are 
carrying out using NMR to probe the bubbles. Neither 
project has reached the stage of detailed agreement with 
first-principles theory, as shown above for the vibrofiu- 
idized bed. Nevertheless some quantitative connections 
with model predictions are made. 



4.1 Measurement of bubble velocity using grain NMR 

Simple and basic questions for the bubbling gas-fiuidized 
bed are: what is the velocity of the rising bubbles, and 
how is the velocity related to the bubble diameter d, gas 
flow rate, and other bed parameters? In analogy with 
gas bubbles in a liquid (or by dimensional analysis), it 
can be argued Vh ~ (dg) 1 ^ 2 , and such a relation has been 
proposed for the gas-fiuidized granular bed with a phe- 
nomenological prefactor [(| . Yet it is not easy to measure 
the mean velocity of bubbles rising in the interior of a 
chaotically bubbling gas-fiuidized bed. 

We have devised a new method to measure the mean 
velocity and other characteristics of the bubbles, using 
NMR on the grains. A pulse sequence (different from 
that shown in Fig. 0) is used to measure the bed density 
profile twice, separated by a time interval At, and the 
experiment is repeated many times. Bubbles in the bed 
appear as localized, downward fluctuations in the bed 
density. When the cross correlation between the fluctu- 
ations in the two time-separated density measurements 
are plotted (Fig.0| a ridge appears corresponding to fluc- 
tuations that propagate at a definite velocity. The mean 
velocity of the bubbles can be accurately measured from 
the slope of the ridge, while the width and shape of the 
ridge presumably indirectly reflect the size and shape of 
the bubbles. This analysis has not yet been completed, 
but it is already clear that the bubble velocity and size 
(from the slope and width of the ridge) are approxi- 
mately consistent with the excess gas flow over minimum 
fluidization flow. Similarly, the bubble size and velocity 
from the graph are approximately consistent with the d 
vs. Vb correlation proposed in Ref. 



4.2 Measurement of bubble-emulsion exchange using 
hyperpolarized gas NMR 

Using the setup shown in Fig. 0] we are able to perform 
NMR experiments directly on the gas within the flu- 
idized bed. Significantly, we have used the same porous 
alumina particles and the same bed diameter as for the 
grain-NMR experiments summarized above. Therefore, 
we can correlate grain and gas NMR information in es- 
sentially identical fluidization conditions (as verified by 
visual measurements of the bed-height hysteresis curve, 
which is extremely sensitive to the fluidization state). 

The 129 Xe NMR frequency is sensitive to the physical 
environment of the gas molecule, which makes it possible 
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Fig. 7 Preliminary measurement of correlations between fluctuations in the density of a gas-fluidized bed of 75-106 /im 
diameter alumina particles, measured using NMR. A small quantity of dodecane is adsorbed onto the porous particles, to 
give NMR sensitivity. The correlations are plotted as a function of the time shift At and the height shift Az. For the left 
plot, there is no gas flow through the bed and hence no fluctuations. For the right plot, there is a gas flow sufficient to put 
the bed into the bubbling regime (superficial gas velocity v g = 1.4 cm/sec). In this case there are large density fluctuations, 
due to the bubbles, and the bubble velocity Vb = 12.4 cm/s can be measured directly from the slope of the ridge in the 
correlation plot. 
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Fig. 8 1 Xe spectra for an alumina-particle bed fluidizcd 
by hyperpolarized xenon gas, measured at 11 different gas 
flow rates ranging from 15 to 125 seem. Three peaks may be 
distinguished: a narrow peak due to gas within the bubbles, a 
broader peak at nearly the same frequency due to gas within 
the interstitial volume of the emulsion (dense granular) phase, 
and a second broad peak shifted to higher frequency due to 
gas adsorbed on the porous alumina particles. As the flow 
rate is increased, the fraction of gas in the bubbles also clearly 
increases. 



to distinguish spectroscopically between gas in the bub- 
bles, gas between the grains, and gas adsorbed onto the 
grains (Fig.|SJ|. Using a suitably tailored NMR pulse se- 
quence |24( , we are able to selectively eliminate the signal 
from gas outside the bubbles and then observe this sig- 
nal re-grow due to gas motion in the bed. In this way an 
important engineering parameter for gas-fluidized beds, 
the bubble-emulsion exchange rate, can be directly and 
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Fig. 9 Preliminary measurement of bubble-emulsion ex- 
change rate as a function of gas flow rate for a bed of alu- 
mina particles fluidized by hyperpolarized xenon gas. The 
exchange rate gives the rate at which gas within the bubbles 
interchanges with gas within the dense granular phase, and 
is a key parameter for engineering applications of fluidized 
beds. The solid curve shows the exchange rate calculated from 
the measured bed expansion along with the phenomenologi- 
cal theory from Ref. Q, while the points with error bars show 
a direct measurement of the exchange rate using 129 Xe NMR. 



noninvasively measured (Fig.[5J). The exchange rate K is 
defined by the following equation [l5| : 

1 dNb Kir r \ 

where Nh is the quantity of a tracer gas within a bubble 
of volume Vb, and Cb,C e are the concentrations of the 
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tracer in the bubble and emulsion phase respectively. For 
our experiments the hypcrpolarized xenon gas serves as 
a completely neutral tracer, and its concentration in the 
various phases can be measured accurately and noninva- 
sively using NMR. 



5 Conclusions 

Using NMR methods, it is possible to measure the den- 
sities and motion of both grains and gas in continuously 
excited granular systems such as the vibrofiuidized and 
gas-fluidized beds. Quantitative results can be obtained 
for many parameters including the grain density pro- 
file, granular temperature (mean random grain kinetic 
energy), grain bubble velocity and size, and gas bubble- 
emulsion exchange rate. 

Where first-principles physical theories exist for the 
granular flow state, e.g., the highly excited vibrofiuidized 
bed, these theories can be directly tested. Conversely, for 
more complex systems like the gas-fluidized bed theoret- 
ical models must make more ad-hoc assumptions (for ex- 
ample that the bubbles behave approximately as bubbles 
in a fluid). In this case non-invasive NMR measurements 
offer possibilities to directly test the results of such as- 
sumptions (e.g., predictions of the bubble velocity-size 
relationship and bubble-emulsion exchange rate). 

In both key goal is to understand the gran- 

ular system sufficiently well to enable extrapolations to 
as yet unexplored regions of parameter space, such as 
reduced-gravity operations. The difficulty and risk of ac- 
tual reduced-gravity experiments render well-supported 
predictions for reduced gravity especially valuable. 
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